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We question the global universe isotropy by probing the alignment of local structures in the cosmic 
microwave background (CMB) radiation. The original method proposed relies on a steerable wavelet 
decomposition of the CMB signal on the sphere. The analysis of the first-year WMAP data identifies 
a mean preferred plane with a normal direction close to the CMB dipole axis, and a mean preferred 
direction in this plane, very close to the ecliptic poles axis. Previous statistical anisotropy results 
are thereby synthesized, but further analyses are still required to establish their origin. 
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INTRODUCTION 

Last years' experiments in cosmology have resulted in 
the definition of a consistent picture of the structure and 
evolution of the universe. The recent data of the cosmic 
microwave background (CMB) radiation, together with 
other cosmological observations, have allowed us to deter- 
mine precise values for the main cosmological parameters 
0. However, the corresponding concordance cosmolog- 
ical model is based on strong hypotheses which need to 
be questioned. They extend from the nature of the grav- 
itational interaction underpinning the cosmological evo- 
lution of the universe, to the physics governing the early 
inflationary era, or also to the cosmological principle for 
the global homogeneity and isotropy of the universe. 

This letter defines an original method to test the global 
universe isotropy through the analysis of the CMB data. 
The observed CMB anisotropies on the celestial sphere 
can be interpreted as a realization of a statistical process 
originating in the inflationary era. The cosmological prin- 
ciple implies the isotropy of the corresponding statistical 
properties. This statistical isotropy has been challenged 
through multiple methods applied to the first year full- 
sky CMB data of the Wilkinson Microwave Anisotropy 
Probe (WMAP) experiment. The detections quoted in 
the following suggest statistical anisotropy with confi- 
dence levels higher than 99 percent. Analyses based on 
N-point correlation functions 0, Q , local curvature Q , 
local power spectra 0, 0, 0, U, and bispectra , sug- 
gest a north-south asymmetry maximized in a coordi- 
nate system with the north pole at (0, cp) = (80°, 57°) 
in Galactic co-latitude 6 and longitude tp, close to the 
north ecliptic pole lying at (9,(p) = (60°, 96°). Analyses 
of multipole vectors, angular momentum dispersion, as 
well as azimuthal phases correlations find an anomalous 
alignment between the low I multipoles of the CMB, sug- 
gesting a preferred direction around (0, ip) = (30°, 260°), 
near the ecliptic plane and close to the axis of the dipole 
lying at (0,<p) = (42°, 264°) [T3iniElElElElEli3. 
Galactic north-south asymmetries are also found in the 



analysis of the kurtosis and the area of the wavelet co- 
efficients of the CMB data. These are mainly due to a 
very cold spot in the southern hemisphere [l^, • First 
results with the angular pair separation method, which 
probes the statistical isotropy both in real and multi- 
pole space, also seem to support those results 20]. On 
the contrary, bipolar power spectra analyses are consis- 
tent with no violation of the statistical isotropy of the 
universe 0, [12, • Finally, theoretical models for an 
anisotropic universe are being studied to account for the 
observed effects |24j . 

Our alternative method probes the statistical isotropy 
of the CMB by the analysis of the alignment of structures 
in the signal. Preferred directions in the universe are de- 
fined as the directions towards which local features of 
the CMB are mostly oriented. The level of preference of 
each direction may be established from simulations, and 
is represented as a signal on the sphere. The approach 
is therefore powerful as it also a priori allows the study 
of the corresponding angular power spectrum in order to 
probe the multipole distribution of the anisotropy. The 
analysis defined relies on a steerable wavelet decomposi- 
tion of the first year full-sky WMAP data. 



DATA AND SIMULATIONS 

The experimental temperature map used for the anal- 
ysis is obtained from the first year WMAP data following 
the procedure originally proposed in [2£( for gaussianity 
tests. First, a best estimation of diffuse galactic fore- 
grounds is removed from each of the eight frequency maps 
of the WMAP receivers at the Q, V, and W bands, as 
prescribed by the specific template fits method in [26q . 
The eight foreground cleaned maps are then combined 
through a noise-weighted linear combination in order to 
enhance the CMB signal-to-noise ratio. The so-called 
KpO mask of |2£| is applied to account for the remain- 
ing strong foreground contamination by diffuse galactic 
emissions around the galactic plane and by bright point 
sources. Finally, the residual monopole and dipole are 
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subtracted outside the mask. All the WMAP data are 
released in the HEALPix pixelization on the sphere [2?| , 
at a resolution identified by the parameter N si de = 512, 
corresponding to maps with several million pixels with a 
spatial resolution of 6.9 arcminutes. The temperature 
map obtained from the above pre-processing is down- 
graded to a lower resolution with N s id e — 32 for the 
specific purpose of our analysis. This provides a map 
with Npi X = 12288 pixels with a spatial resolution of 1.8 
degrees. 

Ten thousand temperature map simulations at the res- 
olution N S ide — 32 are used in order to define confidence 
levels for the analysis results. These simulations specif- 
ically assume the gaussianity and statistical isotropy of 
the CMB. They are generated following the scheme also 
proposed in [23] . A gaussian CMB simulation is obtained 
in spherical harmonics space from the angular power 
spectrum determined by the cosmological parameters of 
the WMAP best-fit cosmological model. The observa- 
tion at each receiver is simulated by convolving that map 
with the corresponding WMAP beam window function. 
After transforming each map to pixel space, a gaussian 
noise realization is also added with the proper dispersion 
per pixel. Following the same procedure as for the data, 
each one of the simulated frequency maps of the eight 
receivers at the Q, V, and W bands are finally combined 
through the noise- weighted linear combination which en- 
hances the CMB signal-to-noise ratio, the KpO mask is 
applied, and the monopole and dipole are subtracted. 

WAVELET FILTERING 

We recall the recently introduced steerable wavelet de- 
composition of a signal on the sphere upon which we base 
our analysis. 

Let the function F(uj) describe a signal on the sphere, 
with the point oj = (0, cp) identified by its co-latitude 
G [0, 7r] and longitude ip G [0, 2ir). A wavelet on the 
sphere is a localized function \I/ that can be dilated at 
any scale a G associated with a given angular opening 
on the sphere, rotated on itself by any angle x £ [0? 27r), 
and finally translated at any position c^o = (#o> <£o)- The 
wavelet filtering of F by \I/ results from the scalar product 
(correlation) of the signal with all the dilated, rotated, 
and translated versions of the wavelet: SP^o^a- Hence, 
the wavelet coefficients W£(luq, X? a ) — (^uj ,x,a\F) char- 
acterize the signal locally at each scale a, in each orien- 
tation Xi an d at each point c^o j3 E3j Individualizing 
the properties of the signal independently at each scale is 
a major advantage of the wavelet scale-space filtering as 
physical phenomena may generally be scale dependent. 
Determining the local orientations of the signal structure 
at each scale also provides essential information for our 
particular analysis. 

Notice that the computation of all wavelet coefficients 



for a large number of scales, and with the same high res- 
olution in local orientation and position is a numerically 
complex calculation. A corresponding fast directional 
correlation algorithm on the sphere was recently defined 
which strongly relies on the use of so-called steerable 
wavelet filters j^. Our analysis of the WMAP data is 
performed independently with the first and second gaus- 
sian derivatives. These steerable filters, represented in 
Mollweide projection in figured are analytically defined 
in 




Figure 1: Mollweide projections of the first (left) and second 
(right) gaussian derivative wavelets at position cjo = (tt/2, 0), 
orientation x — 5 an d scale a = 0.19. Red and blue regions 
respectively correspond to positive and negative values. 



ANALYSIS METHOD 

Here, we consider the analysis leading to the identifi- 
cation of preferred directions from the first year full-sky 
WMAP data at the HEALPix resolution N side = 32. 
The complementary study of the corresponding angular 
power spectrum is postponed to a future work. 

For a given scale a, the angular size of our wavelets 
on the sphere is defined as twice the dispersion of the 
corresponding gaussian. Twelve scales are selected cor- 
responding to angular sizes of the first and second gaus- 
sian derivatives lying between 5 and 30 degrees. At each 
scale, an extended exclusion mask M a is defined in order 
to avoid considering pixels for which the correlation be- 
tween the WMAP signal and the wavelet is contaminated 
by the KpO mask. The pixels added to the initial mask 
are those for which the wavelet coefficients of a constant 
signal in the region of the mask are non-zero, up to a 
given threshold. 

For each pixel uoq of the signal outside the extended ex- 
clusion mask, the direction Xo(^o) for which the wavelet 
coefficient is maximum in absolute value is selected, and 
the corresponding absolute value is retained. This se- 
lects the local wavelet orientation which best matches 
the orientation of the local structure of the signal at each 
point. The great circle is defined which passes through 
that point and admits the corresponding local orienta- 
tion as a tangent. The directions in the sky lying on that 
great circle are considered to be highlighted by the lo- 
cal structure identified, and are weighted by the absolute 
value of the wavelet coefficient |W^(o;o, Xo(^o)? a)\- 

Consider each direction uo on the celestial sphere in 
the set S of the N pix = 12288 pixels at N side = 32. The 
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total weight D a {u) at scale a is the sum of the 7V cros (cc;) 
weights originating from all pixels uJ^ in the original 
signal, with 1 < c < N cros (cu), for which the great circle 
defined crosses the direction considered: 

N cros (uj) 

D «(") = J E \Wi(4 C \xoU c) ),a)\. (1) 

C—l 

The factor A = LN~> E Wo £M a \ W i ( u o, Xo(^o), a) I of 
normalization defines a mean total weight in each direc- 
tion equal to unity for isotropic CMB simulations without 
mask: N~l ^2 UjeS D a (u) = 1. The quantity L = AN side 
stands for the number of points on a great circle on 
a HEALPix grid. Preferred directions in the universe 
are therefore identified as the most weighted ones. No- 
tice that the procedure obviously assigns identical total 
weights to opposite directions. In other words, our direc- 
tions are headless vectors. 

RESULTS 

In this section, the results of the analysis of the first 
year WMAP data proposed above are exposed. The fil- 
tering of the CMB data by the first gaussian derivative 
does not lead to any significant detection. This suggests 
that there is no clear alignment of the local CMB struc- 
tures with the morphology captured by that specific fil- 
ter. The corresponding analysis with the second gaussian 
derivative leads to a strong detection, detailed in the fol- 
lowing. 

At each scale a and at each pixel cj, the total weight 
obtained can be quantified by the number of standard 
deviations a a (uj) through which it deviates from the 
mean /i a (^), as estimated from the simulations. Pre- 
ferred directions, or positive total weights, correspond 
to D a (ou) > [ia(u). Non-preferred directions, or nega- 
tive total weights, are distinguished as D a (uo) < fi a (u). 
We consider in particular the scale as associated with 
an angular size of the filter of 8.3 degrees on the sky. 
The corresponding map of total weights clearly depicts 
the distribution of anisotropy resulting from our analysis 
(figure top panel). 

A rare detection is observed at that specific scale. It 
identifies 20 directions (pairs of opposite points), quali- 
fied as anomalous at 99.99 percent, with an associated 
positive total weight higher than in any of the ten thou- 
sand simulations considered. The 20 anomalous direc- 
tions obtained at scale <23, with positive total weights 
lying in the interval [4.4, 7.4] in a a3 (co) units, are mainly 
concentrated in two clusters around the ecliptic poles 
(figure [3 bottom panel). The most prominent of these 
concentrations comprises 16 of the 20 anomalous direc- 
tions. Its mean position is identified by a northern end at 
(#, if) = (71°, 91°) in galactic co-latitude 6 and longitude 
(p, when each direction is weighted by its corresponding 



total weight. This defines a mean preferred direction in 
the sky with a northern end very close to the north eclip- 
tic pole at (0, (f) = (60°, 96°). The presence of the two 
elongated clusters also suggests that the 20 anomalous 
directions lie on a great circle, defining a preferred plane 
in the sky. The normal direction to the mean preferred 
plane, given by the cross-product of the mean directions 
of each of the two clusters, has a northern end position at 
(#, if) = (34°, 331°). Alternatively, the great circle which 
best fits the 20 anomalous directions identifies exactly the 
same plane. This normal direction lies close to the north- 
ern end of the CMB dipole axis at (0, tp) = (42°, 264°). 

Let us emphasize that 11 directions anomalous at 99.99 
percent are also detected at the neighbour scale cor- 
responding to an angular size of the filter of 10 degrees. 
They belong to two clusters with very similar mean posi- 
tions to those identified at scale as, and therefore single 
out the same structure of anisotropy. 




Figure 2: Mollweide projection of the galactic coordinates 
maps of total weights (top) and directions anomalous at 99.99 
percent (bottom) at the scale ci3 of the second gaussian deriva- 
tive wavelet, and in a a3 units. 



DISCUSSION AND CONCLUSION 

The present original analysis of structures alignment 
in the CMB from the first year WMAP data clearly 
identifies a mean preferred plane in the universe with 
a normal direction close to the CMB dipole axis, and a 
mean preferred direction in this plane, very close to the 
ecliptic poles axis. This result is based on the observa- 
tion of 20 directions anomalous at 99.99 percent. Our 
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original method thus singles out the same directions as 
those highlighted by previous statistical isotropy studies. 
The wavelet approach also identifies the angular size of 
the anomalously aligned structures around 8.3 degrees on 
the celestial sphere, corresponding to a multipole range 
roughly between I = 11 and I = 27. 

This new insight into the anisotropy structure might 
help us to understand its still unclear origin. First, the 
angular size identified for the aligned structures is com- 
patible with the size of primary CMB anisotropics due 
to topological defects such as texture fields [^lf or sec- 
ondary anisotropics due to the Rees-Sciama effect |32j . 
Alignment mechanisms [23, LiJ were recently proposed 
which might be generalized to such structures. Second, 
the identified angular size is also compatible with the 
one above which diffuse galactic emissions dominate the 
CMB emission. Corresponding residual foregrounds can 
be probed by a direct analysis of the WMAP Q, V, and 
W bands frequency maps. Finally, the coincidence of 
the preferred directions detected with the ecliptic poles 
and dipole axes naturally suggests possible unknown sys- 
tematic effects In that respect, notice that the 
angular size of the mesh of the WMAP scan pattern de- 
fined by the combination of the spin and precession of 
the satellite is, again, of the order of several degrees |3fij . 

In conclusion, our analysis provides a first synthesis 
of previous statistical anisotropy results. Further analy- 
ses of the complete total weights distributions at various 
wavelet scales and of the corresponding angular power 
spectra can allow a deeper probe of the anisotropy struc- 
ture. But the origin of the present detection already 
needs to be thoroughly investigated. Various hypotheses 
can be suggested in terms of cosmological or foreground 
structures, or systematics. But nothing at present al- 
lows us to discard the possibility of a global universe 
anisotropy, simple violation of the cosmological principle 
hypothesis. 
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